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Results are presented for simulations of electron-positron 
pair production in relativistic heavy-ion collisions leading to 
electron capture and positron ejection. We apply a two-center 
relativistic continuum distorted-wave model to represent the 
electron/positron dynamics during the collision process. The 
results are compared with experimental cross-section data for 
La 57+ and Au 79+ impact on Gold, Silver and Copper targets. 
The theory is in very good agreement with experiment, veri- 
fying the result that the process increases in importance with 
both collision energy and target atomic number, and improves 
upon previous simulations of this process. 



Early theoretical work on the production of an e~ / e + 
pair through heavy-ion collisions considered only the cre- 
ation of the electron and positron in the continuum. 
However capture by pair production (CPP), in which the 
electron is formed in a bound state of one or other ion, 
becomes a significant process at highly relativistic ener- 
gies. Remarkably this process was sufficiently important 
to enable the synthesis of atomic antihydrogen using the 
low-energy antiproton ring at CERN. A beam of fast an- 
tiprotons impacting on a xenon gas target jjj led to pair 
production with positron capture. Theory had predicted 
HI that cross sections for CPP would increase with en- 
ergy, and indeed this has been verified experimentally 
|§^| . In fact, this process eventually becomes the domi- 
nant mechanism for charge exchange in highly-relativistic 
atomic collisions || . As well as being an interesting area 
of study in its own right, this process has important ap- 
plications in the physics of heavy ion colliders such as the 
large-hadron collider (LHC) and the relativistic heavy- 
ion collider (RHIC) 0. The process of CPP will lead to 
depletion of the charge state of the beam and hence a 
loss in luminosity of the collider. For typical operating 
conditions of such facilities, these losses might amount to 
50% § or more. 

Although the process is strongly coupled at high en- 
ergy, simulations based on relativistic coupled-channel 
calculations |||| have indicated that leading-order per- 
turbation theory is adequate for total cross section esti- 
mates for energies (E) up to 150 GeV/u Q. Nonetheless 
in the energy range E ~ 1 GeV/u where reliable exper- 
imental data exist, theory and experiment have been in 
least agreement. It is this region which we address in this 
paper. 

It is now some thirteen years since Becker and co- 
workers II0N12I obtained the first estimates of cross sec- 



tions for pair creation with simultaneous capture of the 
electron into the ivf-shell of one of the colliding ions. 
However with the exception of Deco and Rivarola, who 
gave a two-center description of the continuum positron 
p3| , two somewhat artificial modes of reaction have been 
distinguished and treated separately when modeling this 
process: excitation from the negative energy continuum 
of an ion to one of its bound states P ,Jll | Jl4| or transfer to 
a bound state of the other ion [ |15|jlfif| . Such approaches, 
while suited to circumstances in which one ion is much 
more highly charged than the other, lack symmetry and 
make a distinction between two separate modes of CPP. 
They lead to different formulae within first-order pertur- 
bation theory [|l5| and hence different projectile charge 
(Zp), target charge (Zt) and E dependencies. As a re- 
sult theoretical estimates of the asymptotic (E —> 00) 
energy dependence of the total cross sections are not in 
agreement, with estimates of ctcpp ~ ln(-E) , and 

more recently |l5| : ctcpp ~ E. The former is based on 
the positron-electron pair being created around the same 
ion, the latter assuming that the pair is divided between 
the two ions. Of course both pathways will interfere and 
contribute to the process, thus pointing to the necessity 
of a two-center treatment for the positron and electron. 
It has been shown |L3| that the two-center description is 
essential in obtaining the correct positron emission spec- 
trum and accurate total cross section for CPP. Leading- 
order perturbation theory (the first Born approximation) 
does give reasonably good estimates for the cross section 
in the high-energy region (E ~ 150 GeV/u) |J for colli- 
sions of heavy ions, and has been a reliable model for fast 
collisions of light ions with low Z targets in the process 
of antihydrogen formation involving CPP by antiprotons 

Experimental results for highly relativistic heavy ions 
on a variety of targets Q support the simple scaling law 
derived from the virtual-photon method (Born approxi- 
mation) which included multiple scattering from the pro- 
jectile ion H alone: ctcpp ~ Z T , for a given energy. At 
lower energies this is not the case |p[ ^8|JTs| ], the Zt de- 
pendence is more complex, showing an enhancement in 
excess of the Z|,-scaling. In this paper we propose a re- 
finement of the Born approximation to take into account 
higher-order scattering processes. In particular we tackle 
the question of the two-center nature of the continuum 
positron and the polarization of the captured electron. 
We find both these effects are vital and lead to theoreti- 
cal results which are in accord with experiment. We dis- 
cuss the physical explanation for scaled cross section en- 
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hancement and provide numerical estimates which agree 
very well with experiment in qualitative and quantitative 
terms. 

Through crossing symmetries the leading-order matrix 
element for the pair production process, in which the 
electron is captured by the projectile P, 



P + T^(P,e~ 



T 



(1) 



is the same as that for the related reaction, e~ +P+T — > 
(P, e~) + T , which is mathematically equivalent to the 
time-reversed ionization process 



(P,e~)+T 



P + T. 



(2) 



In each crossing symmetry the equivalence relies on the 
electron-positron interaction being much weaker than 
their interactions with the highly-charged ions; a reason- 
able assumption. Let rp,t and r' T ,t' be the space and 
time coordinates of the electron in the projectile and tar- 
get frames, respectively. The nuclei follow straight-line 
paths with relative velocity v. The Hamiltonian, in the 
projectile frame of reference and in atomic units, is given 
by: 

H = -ica.V rp + f3c 2 + V P (r P ) + S 2 V T (r' T ) (3) 

where a and (3 are Dirac matrices and S is the operator 
which transforms the wavefunction from the projectile 
frame to the target frame, namely 



S= (i + 37)^(1 - xct.v), 



(4) 



with x = U7c _1 (7 + 1) -1 ,7 = (1 — u 2 /c 2 ) -1 / 2 , and 1 
represents the unit matrix. For a given impact parameter 
b, the transition amplitude can be written in the form EG] 



A(b) 



dt I dri 



x\{H-id t )Xi 



(5) 



where x% ancl Xf are the initial and final states. 

The undistorted bound-state is approximated by a 
semirelativistic {Zx <C c) wavefunction : 



Oi 



where 



and 



$ H = (2ic) 1 a.V rp $oi, 



(6) 



(7) 



(8) 



with E s i the non-relativistic eigenenergy, and the elec- 
tron spin along the beam axis defined as 'up' by: ujf — 
(1 0) and 'down' by uj = (0 1 0) 
The continuum function is given by: 



o/ 



(9) 



where 

$0/ = {2tt)-^ N*{uj P ) 1 F 1 {-iujp;l; -ij e (v e r P + v e ■ r P )) 

The spinor correction term is given by 

$1/ = {2TT)~i{2i leC y 1 N*{u J p) x 

a.V rp iFx(— iuj P ; 1; -i7 e (u e r P + v e .r P )) 

x e-^t+^^Svluif. (11) 

with top — Zp/v ei where v e is the electron velocity. 
N(Q = exp «/2) T(l - iQ and, 

S Ve = (5 + \le)%{l-x e a.v e ), (12) 

where a; e = VejeC 1 ^ + l)" 1 , and 7 e = (1 - d 2 /c 2 )~ 1/2 . 
These functions are appropriate when Zp_x <S c. 
The initial distortion factor L i is a matrix given by: 



where 



and 



Li — L 0i + L u , 



L' 0l = exp(— ivr\n.[yvr T + -fv ■ r' T ])l, 



L' u = S- 1 (2i 7 c)- 1 a.V rir L 0i S, 



with vt = Zt/v. 

The final state distortion is given by |TJ 

L f = L 0f + L lf> 



(13) 
(14) 
(15) 

(16) 



where 

L' of = N*(ii}' T )iFi(—i(j T ; 1; -ij' e (v' e r' T + v' e ■ r' T ))l (17) 
and 

L' lf = S- 1 (2 ll ' e cr 1 a .V r , T L' of S. (18) 

Retaining terms of first order in Z/c, we have relativis- 
tic continuum distorted wave eikonal initial state (RCD- 
WEIS) wavefunctions HQ: 

Xl = L' 0l $ 0l + L' u <$> 0i + L' 0i $ u , (19) 



Xf = L' of $ of + L' lf $ f + L' of $ lf . 



(20) 
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FIG. 1. Scaled cross sections, ocppIZ t in microbarns, for 
pair-production with electron capture by fully stripped Lan- 
thanum ions (La 57+ ) striking thin foils of Copper (Z T = 29), 
Silver (Zt = 47) and Gold (Z T = 79). Comparison with 
RDWB theory for capture to the ls-state. 
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The triply differential cross section, with respect to the 
electron momentum (p e ) is defined as: 

a( Pe ) = (da CPP /d Pe ) = J db\A(b)\ 2 . (21) 

Using the Fourier transform method pfj] we define 



T(rj) = -fv J db exp(-iri ■ b)A(b), 



(22) 



where T(rj) is a product of single-center integrals. The 
total cross section is obtained from the integral over the 
ejectile momentum (or velocity) and takes the form: 



ccpp 



E 

spins 



i 



1 1 d0sm9 
o 



d V \T( V )\ 2 . 

(23) 



where we sum over all the spin states of the electron and 
positron pair. 

In order to compute CPP cross sections (|]) we note 
that a positron with energy e + and momentum p + trav- 
eling forward in time in the final state is equivalent to 
an electron with energy — e + and momentum — p+ in the 
initial state. Thus we must take: 



(qrlL/z/o 

FIG. 2. Scaled cross sections, acpp/Z T in microbarns, for 
pair-production with electron capture by fully stripped Lan- 
thanum ions (La 57+ ) striking thin foils of Copper (Z T = 29), 
Silver (Z T = 47) and Gold (Zt = 79). Comparison with 
RCDWEIS theory for capture to the ls-state. 

We compare our results for the Relativistic First Born 
(RIB), projectile centered, approximation, in which the 
initial and final state distortions are neglected, and the 
Relativistic Distorted Wave Born (RDWB) approxima- 
tion [ jl3| , where the two-center positron wavefunction is 
used but the initial state distortion is omitted. The Born 
approximation, which assumes that the positron is in the 
continuum of only one of the ions, is ambiguous. The 
Born approximation of Bertulani and Baur || assumes 
the positron to be a projectile-centered wavefunction, 
while the Born approximation of Eichler Jl5| takes the 
positron to be centered around the target. These two 
models can be viewed as approximations to the wave- 
function ( p0| ) in which lu' t = and ujp = 0, respectively. 
By retaining both scattering center contributions the in- 
terference effects are taken into account. In comparing 
RDWB and RIB, it is known that these two-center in- 
terference effects reduce the cross section for CPP in the 
relativistic domain jL3|. This suppression of CPP is the 
converse of the two-center enhancement (capture to the 
continuum) that arises in ion-atom ionization pC[ | and is 
analogous to the effect of the Fermi function for (3^ decay 
ill. 



v e 



v e 



(24) 



The experiments of Bclkacem et al. ||,|T^,|l9| were 
for fully stripped Lanthanum ions (La 57+ ) striking thin 
foils of Copper (Z T = 29), Silver (Z T = 47) and 
Gold (Zt = 79). The collision energies were E = 
0.405,0.956 and 1.300 GeV/u. The two graphs pre- 
sented compare the scaled total cross sections (ctcpp /Z~t) 
given by theory and experiment. Consider figure |l| which 
compares RIB and RDWB with the measurements. Of 
course, the scaled RIB curve is independent of Zt, and it 
clearly shows the increase in importance of CPP with in- 
creasing collision energy. Considering the RDWB model 
however, we see a progressive reduction in the scaled cross 
section as Zt increases. This is in agreement with the 
findings of Deco and Rivarola [^3|, who reported a de- 
crease in the size of the singly differential cross sections 
by an order of magnitude. Their model is essentially 
our RDWB approximation. While this model shows Zt 
dependence for the scaled cross section, the trends and 
absolute values are incorrect. It predicts a suppression of 
the scaled cross section rather than an enhancement as 
Zt increases. Thus the RDWB theory data for Gold 
gives the lowest scaled cross-section while experiment 
shows that it should be the highest. The same incorrect 
trend was obtained by the target-centered Born approx- 
imation |l5f| . 

In contrast (figure ^) the equivalent results for RCD- 
WEIS show the observed enhancement with increasing 
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Zt- However the theoretical data lie below the experi- 
ment for the more energetic collisions. In comparing with 
experiment we have only presented simulations for the 
dominant channel; capture to the Is ground state. At 
very high energies capture to excited states is thought 
to contribute ~ 30% to the total capture cross section 
]6|.p2[ . This would partly explain the differences between 
our theory results and the experimental data. These re- 
sults are a great improvement on the RIB approximation 
and are consistent with the experimental data. 

Other experimental results are available for the impact 
of faster and more highly charged beams: 10.8 GeV/u 
Au 79+ @ and 0.956 GeV/u U 92+ (§ for the same tar- 
gets. The Gold beam results (table [i]) indicate that the 
Zf. dependence is established at the higher energies, as 
predicted by the simple projectile-centered Born approx- 
imation [ |10| . Even at this higher energy our theoretical 
results (table ||) show an enhancement in excess of Zj.. 
The experiment is in much better accord with the flat 
scaled cross section data given by the Born approxima- 
tion g(J. 

Finally we note that for very large charges, for ex- 
ample U 92+ , the semirelativistic approximations for the 
wavefunctions are not valid. Qualitatively this can lead 
to underestimation of the cross sections for high charges 
when using approximate wavefunctions [g_ljjl5j|. In such 
cases, Coulomb-Dirac wavefunctions are required in or- 
der to model the process [ [Tl|| . 

In summary, we have proposed and tested a new dis- 
torted wave model which incorporates and reduces to ap- 
proximations used previously to describe CPP. We con- 
firm that, as has previously been shown [fl3"[ , the inclusion 
of distortions from both ions on the positron continuum 
state leads to a reduction in the cross sections. However 
including distortion of the bound electron leads to an 
increase in the total cross sections and a more accurate 
fit to the experimental data for fully-stripped relativis- 
tic Lanthanum ions. This improvement in results when 
including distortions demonstrates once more the neces- 
sity of a two-center treatment for an accurate theoretical 
description of this reaction. However our cross section 
predictions for faster and more highly charged Gold ions 
do not accord with the experimental data which show a 
Z\ dependence. While the refinements introduced in our 
model are significant theoretical improvements, clearly 
there still exist several unresolved important differences 
between theory and experiment. 
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TABLE I. Total cross sections, <tcpp in barns, for electron 
capture from pair production for 10.8 GeV/nucleon Au 79+ 
impact on Gold, Silver and Copper foils. 



Zt 


Experiment [7 


CDWEIS theory 


Becker et al 1C] 


79 


8.8 ± 1.5 


15.85 


10.1 


47 


4.4 ± 0.73 


3.44 


3.6 


29 


1.77 ± 0.31 


0.74 


1.36 



R.J.S. Lee and J.V. Mullan acknowledge financial sup- 
port from the Department of Education for Northern Ire- 
land through the Distinction Award scheme. 



4 



